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We study neutrino masses and mixing in Supersymmetric Models without R parity and 
with generic soft Supersymmetry breaking terms. Neutrinos acquire mass from various 
sources: tree level neutrino-neutralino mixing, loop effects and non-renormalizable oper- 
ators. Abelian horizontal symmetries (invoked to explain the smallness and hierarchy in 
quark parameters) replace R parity in suppressing neutrino masses. We find lower bounds 
on the mixing angles: sin6'ij > m{£~)/m{£~) (z < j) and unusual order of magnitude pre- 
dictions for neutrino mass ratios: m{i>f,) /m^u^) ~ sin^6'i2; m{i>i) / m{v-j-) ~ 10~^sin^6'i3 
[i = 1,2). Bounds from laboratory experiments exclude rn^j^ > 3 MeV and cosmological 
constraints exclude m,^^ > 100 eV. Neither the solar nor the atmospheric neutrino prob- 
lems are likely to be solved by — Ue oscillations. These conclusions can be evaded if 
holomorphy plays an important role in the lepton Yukawa couplings. 
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1. Introduction 



The search for neutrino masses is one of the most promising directions to find evidence 
for the incompleteness of the Standard Model. Theoretical input is required in order to 
direct experiments to the most plausible values of neutrino masses and mixing angles. In 
particular, an understanding of the neutrino sector by the same means that explain the 
quark and charged lepton parameters would be desirable. Supersymmetry combined with 
horizontal symmetries can provide this understanding 

In Supersymmetric models with the MSSM particle content and with R parity (Rp), 
lepton number is violated by non-renormalizable terms only. Terms of the form j^L(j)uL(j)u 
(M is a high energy scale, L is the lepton doublet and (pu is the Y = +1/2 Higgs doublet) 
lead to neutrino masses of the see-saw type, rui, ~ The consequences of Abelian 

horizontal symmetries in this framework were investigated in ref. [§. A number of inter- 
esting order of magnitude relations among the lepton parameters were found to hold in a 
large class of models: 

^~sin2%, (1.1) 
sin».,>^, (1.2) 

^>(^)\ (1.3) 

where i < j and Ug^^u^^UT denote the mass eigenstates with mixing of 0{1) with e, /U, r, 
respectively. Interestingly, predictions analogous to (|1.2|) and ( pT^) apply to the quark 
sector (namely, Vij > and Vckm ~ 1) and are experimentally valid g. 

The same horizontal symmetries that explain the smallness and hierarchy in fermion 
parameters can naturally solve the problems related to lepton flavor and lepton number 
violation that arise in Supersymmetric models without Rp [0-PJ^.0 In this case, lepton 



While horizontal symmetries can rather easily take the role of Rp in suppressing lepton 



number violation, it is much more difficult to do so for baryon number violation |1C]. Therefore, 
as in ref. Q , we simply assume that baryon number is a symmetry of Nature. 

1 



number is violated by renormalizable terms, and the resulting phenomenology is strikingly 
different from the one predicted by Supersymmetric i?p-symmetric models. In this work 
we examine the question of lepton masses and mixing angles in models of Supersymmetry 
without Rp but with a horizontal symmetry. 



2. The Theoretical Framework 

We work in the framework of the Abelian horizontal symmetry Ti. that has been 
introduced in refs. ||Tl| - p!3| . Ti is explicitly broken by a small parameter A to which we 
attribute charge -1 (and a numerical value of C(0.2), to explain the Cabibbo angle). 
This can be viewed as the effective low energy theory that comes from a Supersymmetric 



extension of the Froggatt-Nielsen mechanism at a high scale [^. Then, the following 
selection rules apply: 

(a) Terms in the superpotential that carry charge n > under TC are suppressed by 0(X'^), 
while those with n < are forbidden by holomorphy; 

(b) Terms in the Kahler potential that carry charge n under Ti. are suppressed by 0(AI"I). 
Without Rp (or lepton number), there is a-priori no distinction between the Y = —1/2 

Higgs doublet (pd and the three lepton doublets Lj. (Wherever convenient, we denote the 
four doublets by Lq,, ct = 0, 1, 2, 3.) The four doublets, however, carry in general different 
horizontal charges H. We identify the Higgs doublet with the doublet field that carries the 
smallest (positive) charge, which we choose to be Lq (we use interchangeably Lq = ^d), 
and we order the remaining doublets according to their charges: 

H{Li) > H{L2) > H{Ls) > H{<Pd) > . (2.1) 

A similar ordering is made for the three generations of ii (charged lepton singlets), Qi 
(quark doublets), and di (down quark singlets). The model is phenomenologically viable 
if (for tan/3 ~ 1) the following condition holds [§]: 

H{L,)-H{<Pd)>3. (2.2) 

Our methods of analyzing lepton and neutralino mass matrices are described in detail 
in refs. H and M , respectively. Specifically, we use the following selection rules to estimate 



the magnitude of the various contributions. For the quadratic terms in the superpotential, 
lia.Loc4>u, it is : 

^ \ mX\^iL.)+Hick^)\ H{L^) + H{(t)u) < 0. ^ > 

Here jl is the natural scale for the terms and m is the Supersymmetry breaking scale. For 
simplicity we assume that fl is 0{m). As /xq is phenomenologically required to be of 0(m), 
we take H{(pd) + H{(pu) ~ 0. Modifications to the case where the natural scale for /i is, 
say, Mpianck and it is suppressed down to rh by the horizontal symmetry, as in the models 



of [13| and Q, are straightforward. The selection rule for the coupling of the quadratic 
soft Supersymmetry breaking terms, BaLct(j)u (here stand for the scalar components) 
is: 

S„~m2Al^(^'^)+^('^")l. (2.4) 
Finally, the selection rules for the trilinear terms X'^-j^L^Qjcik and XapkLaLplk are 

y f XHiL^)+H(QMH{d.) H{L^) + H{Q,) + H{dk) > 0, ^1 
^ I H{L^) + i7(Q,) + H{dk) < 0, ^ • ^ 

A ~ I A^(^")+^(^/')+^(^"'=) H{L^) + if(L^) + if (4) > 0, , . 

10 H{L^)+H{Lp)+H{h)<^. 
Note that Aq^^^ and Xojk are practically the Yukawa couplings for the down sector and for 

the charged lepton sector. 

The order of magnitude relations ( p..l| )-( P~3[ ) were derived in a large class of models 

where all entries in the lepton mass matrices carry positive charges. (They are actually 

applicable in a larger class of models, where the holomorphy-induced zero entries, if any, 

do not affect the physical parameters.) In this work we restrict ourselves to this class of 

models. 

3. Neutrino Masses and Mixing 

There are several important sources for neutrino masses in this framework, each giving 



a different scale: renormalizable tree-level mixing with neutralinos [|T5| - P4[ ; quark-squark 



and lepton-slepton loop corrections ||T6| , p5l - p2| ; and non-renormalizable see-saw contribu- 



tions [P3| - p4| . We now discuss each contribution in turn. In our various estimates we take 
tan /3 ~ 1. 



(i) Renormalizable tree-level contributions. 
These contributions arise when the ^-terms in the superpotential and the Supersymmetry 
breaking B terms in the scalar potential are misaligned, Ba 7^ -B/Ua , or the Supersymmetry 
breaking scalar masses do not satisfy the eigenvalue condition ml^pup = fh? H- This 
yields misalignment between the VEVs Va = {L^) and the terms, 

Sin 4 = , (/i = HaHa , = Vo,Vc) (3.1) 

which induces neutrino mixing with the neutralinos. Only one neutrino acquires a mass 
from this effect: ~ mz sin^ ^ (here mz stands for the electroweak or Supersymmetry 
breaking scale). In the absence of any symmetry reason for alignment, we expect sin^ ~ 1 
and the natural scale for m,^ is the electroweak scale. The further required suppression 
comes from 7Y-violation, Vi/vo ~ Xl^(^i)-^('l'd)] g. Thus, (|3TD together with (|2TD gives 



The massive neutrino is then Uj- , which is close to the interaction eigenstate with the 
smallest horizontal charge among the Li. The experimental upper bound m,^^ < 24 MeV 
||35|| , when confronted with ( |3.2[ ), is the source of the constraint ( p.2| ). 
(ii) Quark-squark loop contributions. 
Loops with down quark and squarks contribute 

1 nju loop '^Kkl'^'jmn {M'^)lm{Mlj^)kn 

2 Idtt^ 

where M*^ is the d-quark mass matrix, and Mf^j^ is the left-right sector in the (i-squark 
mass-squared matrix. The experimental value of Vcb ~ ms/rrib strongly suggests that 
H{d2) = H{d3) and consequently (M<^)32 ~ (M^)33 ~ and (Mf^)32 ~ (^^2^)33 ~ 
mnib. From this, together with ( |2.5|) , we learn that the largest contributions to ( |3.3| ) come 
from k = m = 3, and / , n = 2, 3. This, in general, gives mass to the two light neutrinos. 
Taking into account that A033 ~ rrih/mz, one obtains: 



m 



loop 



.loop 



>op>2[H(L2)-H(L3)] 



^v.- . . ^loop^^2[H(Li)-H(L3)]^ (3.4) 

1 ?>rnt _ T 



which imphes the upper hmit for the loop-induced masses m,^- < IQ-^m^^ <leV,i = 1,2. 

Note that if the Supersymmetry breaking trihnear scalar couplings are proportional 
to the Yukawa couplings {Mlr ~ wM*^), the dominant quark-squark loop contributions 
to (|3.3| ) yield a degeneracy in the mass matrix, and only one neutrino in ( |3.4| ) acquires 
mass. In the presence of a tree-level mass ( p.2|) for Ur , this mass eigenstate is z^^ (and 
is close to the interaction eigenstate L2). The same result applies also to the case that 
(M'^)32 < (M^)33 and {M 1^)32 < {Mfj^)33. Then a contribution to mjy^ arises from 
quark-squark loops with e.g. /c = n = 2, / = m = 3, giving m^^/m,^^ ~ ^^m^ ^2[h(Li)-h(L3)]^ 
about two to three orders of magnitude below ( |3.4| ). This is somewhat smaller than the 
contribution from lepton-slepton loops discussed below, 
(iii) Lepton-slepton loop contributions. 
Loops with charged leptons and sleptons contribute 



1 1,^1' loop Kkl^jmn {M'^)lm{M^j^)kn 



2 16^^2 • ^ • I 

Using (|2.6| ), we learn that the largest contribution from ( |3.5D has 



,4 



= - 10-^ (3.6) 



about two orders of magnitude lower than the dominant quark-squark contributions. As 

already mentioned, this contribution plays a significant role only when Mf^^ — rhM'^ holds 

to a good approximation. 

(iv) Non-renormalizable contributions. 

The dimension-5 terms -^i'ii'j4>u4>u give 0: 



2 M 

This, in general, contributes to both light neutrinos: 

^ ~ gnr^2[//(L2)-H(L3)] 



(3.7) 



n T 



m 
m 



nr 



!^ ^^^X^[HiLi)-H{L3)]^ (3.8) 



M V M 



The relative importance of the non-renormahzable and loop contributions to m^^ and 
rrii,^ depends on the scale M (which is, roughly speaking, the natural scale for the masses 
of right-handed neutrinos). For M > 10^ GeV ^ the leading contributions come from loops, 
while for M < 10^ GeV, the non-renormalizable contributions dominate. 

Adding up the various contributions, and defining 

e = max(e"^e^°°P), (3.9) 

leads to the following order of magnitude estimates for the neutrino masses and mixing 
angles: 

/mz~ A2[^(^^)-^(<^^)], 

(3.10) 

m.7m., ~ eA2[^(^^)-^(^3)] (i = l,2), 

sin% ~ XHiL.)-HiLj) ^ ^3_;^^^ 

The charged lepton mass ratios are estimated to be 

m,Jm,^ ~ ;,H(L.)+H(£-0-H(L,)-H(/,) _ ^3_;L2) 

The charged current mixing matrix mixes not only the leptons among themselves, but 
also leptons with higgsinos and gauginos: 

(3.13) 

sin0,^^^~A^(^^)+^(^-). 

The fact that the neutrino mass eigenstates have an isotriplet ws component in them, leads 
to flavor changing couplings of the Z-boson to neutrinos ~ gQijZuiUj-. 

The estimates ( |3.1C1| ) and ( |3.11| ) give the following relations between neutrino masses 
and mixing angles: 



nil, 9 
— - ~ esm t/23, 

. 2 



esin ^13, (3.15) 



^ sin^ 6*12 
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There are two order of magnitude relations that are independent of e : 

~ sin6>i2 ~ ■ ^ , (3.16) 

rrij,^ sm 6*23 

so that for the two hght neutrinos ( |1 . 1| ) stiU holds. The order of magnitude inequality ( |1.2| ) 
is maintained: 

sin^,, >^ (z<j). (3.17) 
The relation (|1.4| ) is also maintained, 

rriu, <mu^ <m^^. (3.18) 



However, unlike (|1.3| ), the light neutrinos are much lighter than v^: 



m 



^ < e (^ = 1,2). (3.19) 



For a scale M > 10^ GeF, this gives — — < 10 ^. It is interesting that these models can 



naturally give mixing angles of (9(1) with the third generation |2^] while the corresponding 



mass ratios are very small. (For different mechanisms that give such a result, see [ PB| -| 



4. Theory Confronts Experiment 

In this and the next sections we show that the order of magnitude relations derived in 
the last section, when combined with various experimental and cosmological constraints, 
exclude large regions of the mass-mixing parameter space. Most of our discussion in these 
two sections is independent of the question of Rp violation. 

As the charged lepton masses are known, eq. ( |3.17| ) provides significant lower bounds 
on the lepton mixing angles. With rrie/rn^ ~ and m^/rrir ~ A^, we get 

sin^23>A2, sin0l3>A^ sin^i2 > A^. (4.1) 

The lower bound on sin ^23 is particularly significant. First, if m,^^ is in the appropriate 
range, z/^ — Uj- oscillations will be observed in the CHORUS, NOMAD and E803 experi- 
ments. Second, combining it with the upper bound sin^ ^23 = BR{tx fii^r) < 6 x 10~^ 



for m,y^ > 3 MeV ||3^ results in 

m,^ < 3 MeV. (4.2) 



Third, in combination with the bound on — Vj- osciUations, sin^ 2^23 ^ 0.004 for 
Am^ > 100 it givesi 



sin ^23 ~ for 10 eV < m^^ < 3 MeV. (4.3) 

As we predict sin 6*1 3 < sin 6*23, (^3) imphes also 

sin^i3<A2 for 10 < m,^ < 3 MeV. (4.4) 

This bound is stronger than the bound from z/e — Vt oscillations, sin^ 2^13 < 0.12 for 
Am^ > 100 eV'^ ||39|| , which, in this range, gives sin0i3 < A. The latter bound, however, 
holds independently of whether holomorphy plays a role in determining the mixing angles. 

Eqs. (f4.2|), ( ^7^ ) and ([4.4|) are applicable also in models with because they result 
from ( |3.17| ) which holds independently of violation. 

5. Theory Confronts Cosmology 

Cosmological considerations related to the age and the present energy density of the 
Universe provide a constraint on the mass and lifetime of neutrinos. For masses in the 
range 100 eV - a few MeV, the constraint reads (see e.g. |l4T[| ) 

ml^r^^ < 2 X 10^° eV^ sec. (5.1) 

The framework of Abelian horizontal symmetries allows an estimate of the neutrino decay 
rates. It is interesting to find whether could have a fast enough decay mode to fulfill 
( p7l| ) and have its mass above 100 eV. 

The dominant decay modes are most likely those which proceed via gauge interactions. 



The bound (|4.2| ) leaves only a very small window where the VF-mediated tree level z/^ 



e+e Uf, is allowed. The rate can be estimated to be: 

r(z/, ^ e+e-z/e) ^t.2n 

— sm 6/13 



^ ' 2 (5.2) 



^1 r.^ ~ 3 X 10^1 eV^ sec. 



^ The interpretation of oscillation experiments might change in the framework of Supersym- 
metry without Rp because new neutrino interactions are introduced In our case, however, 

the horizontal suppression makes these new interactions practicaUy negligible. 



Together with ( [4 .21) , we find that (|5.1|) is satisfied only for sin ^13 > A^. Since there are 
charged particles in the final state, however, a stronger bound (from considerations of the 
cosmic microwave background radiation) applies, Tj^^ < 10^ sec. This cannot be satisfied for 
m^^ < 3 MeV and sin 6*13 < A^. Furthermore, detailed studies of the effects of a massive i^r 
during the Nucleosynthesis era [^2| - |4^ suggest that for m^^ > 0.5 MeV, and independently 



of the decay modes, t^^ < 10^ sec is required, which closes the window even more firmly. 
Therefore we conclude that Ut- —>■ e~^e~Ue does not open any window for a heavy Ut- Again, 
this conclusion holds also for models with Rp. 

All other decay modes are fiavor changing neutral current processes. There are three 
types of contributions to such processes: 

(a) Loop diagrams with gauge particles, suppressed by the charged current mixing angles; 

(b) Tree level Z-mediated decays suppressed by the O^j mixing angles; 

(c) Tree level slepton-mediated decays suppressed by the selection rules for the Ajj^ cou- 
plings. 

The first class is common to models with and without Rp, but the other two are 
present only in i2p-violating models. In any case, we found that none of these channels is 
fast enough to allow m^^ > 100 eV. For example, the rate for the Z-mediated Uj- 3z^^ 
can be estimated to be: 

^ ' / 5.3 

10^3 



023 

This is significantly suppressed compared to ( ^.2|) and does not satisfy (^T|). 

As Vj- is predicted to be the heaviest among the neutrinos, we conclude that in the 
framework of Supersymmetry and Abelian horizontal symmetry with or without Rp (and 
assuming that holomorphy does not play a role in determining sin ^23) 

m,, < 100 eV (5.4) 

holds for all neutrino masses. We note that in the framework of a single or broken 
by A ~ 0.2, this requires H{Li) — H{(pd) ^ 6 which may be too large for reasonable models. 
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In some models of ref. |11], however, where the symmetry breaking parameters are much 
smaUer, this can be achieved with charge differences < 2. 



6. Solar and Atmospheric Neutrinos 

The upper bound m^^ < 100 eV leads to even stronger bounds on m^^^ and m^^ . These 
bounds, however, depend on M and on tan (3. We believe that the most likely situation 
is (a) M > 10^ GeV (which, for example, applies in all models where the Supersymmetric 
extension of the Standard Model is valid up to some GUT scale) and (b) tan /3 ~ 1 (which 
is the natural value [^5]-^). Then in ( |3^ ) e ~ 10~^ leads to 



m,^ < 100 eV, m^^ < m,^ < 10"^ eV. (6.1) 

This bears important consequences for the solar and atmospheric neutrino problems. The 
value m^^ < 10~^ eV is inconsistent with 

ml^ -ml^r^Qx IQ-^ eV^ (6.2) 

that is required to solve the solar neutrino problem through the MSW mechanism (see e.g. 
Il47[]), but is consistent with m^^^ — m^^ ~ 10~^^ eV^ that could solve it through vacuum 
oscillations (see e.g. [0). It is also inconsistent with 

ml^ - ml^ ~ 10-2 eV^ (6.3) 

that is required to solve the atmospheric neutrino problem through — Uf, oscillations 
(see e.g. [p| ). 



The solar neutrino problem can still be solved by the MSW mechanism with Uj- — v^, 
oscillations (if m^^ ~ 10"^ eV^ or the atmospheric neutrino problem can be explained by 
" oscillations (if m^^ ~ 10"^ eV). However, we find that: 

a. The two problems cannot be solved simultaneously; 

b. The required horizontal charges are inconveniently large: H[L^) — H{(pd) ~ 10 and 
H{Li) — H{(pd) ~ 12 to get the small angle MSW solution for the solar neutrino 
problem, and H{L2,3) — H{(f)d) ~ 8 to solve the atmospheric neutrino problem; 
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c. Such a light Uj- does not contribute to the dark matter and cannot play any role in 
structure formation. This would require 



m^^ ~ 10 eV. (6.4) 



This situation is very different from the Supersymmetric models with i?p, where 
and ( |6?^ ) can be simultaneously accommodated [Q. 

Things are different if we relax either of our two extra assumptions. As the tree 
level contribution to m^^ is suppressed by tan^ (3 and the loop contributions are enhanced 
by tan"^/9, a large tan/3 would give a large e {e.g. tan/3 ~ 20 gives e ~ 10-2). Then 
we can easily accommodate the dark matter ( |6.4| ) and solar neutrino ( |6.2| ) constraints 
[H{L2) = H{L^) + 1 = H{(pd) + 6). Alternatively, the solar and atmospheric neutrino 
problems can be solved simultaneously {H{Ls) = H{L2) = H{(f)d) +6). It is non-trivial 
that, in this scenario, H{L2) and H{L^) which are fixed by the requirements on m^^ and 
rrijj^ give, at the same time, sin ^23 = C^(l) cis required to solve the atmospheric neutrino 
problem. 

For M < 10^ GeV, the non-renormalizable contributions to dominate over the 
loop corrections. In this case m^^ can account for the dark matter, while m^^ can accom- 
modate the solar neutrino constraint ( |6.2| ) (however, this requires M < 10^ GeV). 

Finally, if holomorphy does play an important role in the physical parameters, then 
even the conclusions of sections 3-5 can be evaded. For example, we can construct models 
where sin ^23 ^ A^, which would allow for m^^ above the 3 MeV bound of section 4. Then, 
with sin ~ (which is marginally compatible with the experimental limits [pO| - |5T| ), 



the decay Uj- e~^e~Ue can still open a window for m^^ close to its experimental bound. 
Explicit examples of models where holomorphy induces approximate zeros in the mass 
matrices and affects physical parameters can be found in refs. [jT2|,|T^,|[. 



7. Discussion 

Models of Supersymmetry with Abelian horizontal symmetries have interesting impli- 
cations for neutrino masses and mixing. We distinguish three cases: 
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(i) Models with Rp. 

(ii) Models without Rp and with generic soft Supersymmetry breaking terms (up to the 
selection rules from the horizontal symmetry) . 

(iii) Models without Rp but with universal Supersymmetry breaking terms implying align- 
ment at a high scale (namely Sa = B/ia and m^^ = fh?5a(3)- 

Class (i) was analyzed in 0. Class (ii) has been studied in this work. Class (iii), 
which yields a scenario quite different from the one investigated here, will be discussed in 
a forthcoming paper. We now compare the predictions of class (ii) with those of (i). 

In a large class of models, where holomorphy does not introduce zeros in the mass 
matrices (or, if there are such zeros, they do not affect the order of magnitude of the 
physical parameters), we find the following order of magnitude relations among the mass 
ratios and mixing angles: 

m^Jm„^ m^Jm„^ 

(z) sin^^as sin^^is sin^ 6*12 (7.1) 

(ii) 10-^sin2^23 IQ-^sin^^ig sin^ ^12 

Note that the following relation among the mixing angles holds in both classes: 

(z) , (n) : sin ~ sin 6*12 sin 6'23- (7.2) 

Furthermore, since sin^^j > mi. /mi. holds independently of Rp, 

(z) , [ii) : sin^ 023 > 10-^ sin^ > 10"^ sin^ 612 > 10"". (7.3) 

This leads to the following ranges for the mass ratios: 

m^^./m^^ rui^jm^^ m,yjm^^ 

(i) 10-3 - 1 10"^ - 1 10"^ - 1 (7.4) 

(ii) 10-1° - 10"^ 10"^^ - 10"^ 10"'^ - 1 

We conclude that measurements of the lepton mixing angles would test the Supersymmetric 
Abelian horizontal symmetry framework while measurements of neutrino mass ratios will 
serve to distinguish between models with or without Rp. 
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In ref. it was shown that m^^jmjj^ > {m^/m-r)'^ together with cosmological 

considerations, strongly suggests that all neutrinos are lighter than (9(100 eV). Models 
without Rp predict m,^^/m^^ ^ {m^/rrir)'^ but we still find that all neutrinos are lighter 
than (9(100 eV). This is a consequence of the fact that there is no decay mode large 
enough to fulfill the cosmological constraints on massive neutrinos. 

In models with Rp, one can accommodate m^^ to contribute sizeably to the cosmo- 
logical dark matter, as well as m,j^ in the correct range required by the MSW solution of 
the solar neutrino problem. In models without Rp (and without any alignment condition) 
we find that, unless the scale of New Physics M is surprisingly low ( < 10^ GeV), is 
too light to play any role for matter enhanced oscillations of the solar UeS. 

Finally, we emphasize that our various predictions are not entirely generic to models 
of Abelian horizontal symmetries. As described briefiy in section 6, a large tan/3 and/or a 
small scale M would modify our discussion of the solar (and atmospheric) neutrino prob- 
lem. But more important, one can construct models where holomorphy plays an important 
role and circumvents the otherwise model-independent predictions of eqs. (|3.16|) , ( |3.17| ) 
and ( CT ). 
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